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Landau damping of Gardner solitons in a dusty bi-ion plasma 
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The effects of linear Landau damping on the nonlinear propagation of dust-acoustic solitary waves 
(DASWs) are studied in a collisionless unmagnetized dusty plasma with two species of positive ions. 

The extremely massive, micron-seized, cold and negatively charged dust particles are described by 
fluid equations, whereas the two species of positive ions, namely the cold (heavy) and hot (light) 
ions are described by the kinetic Vlasov equations. Following Ott and Sudan [Phys. Fluids 12, 2388 
(1969)], and by considering lower and higher-order perturbations, the evolution of DASWs with 
Landau damping is shown to be governed by Korteweg-de Vries (KdV), modified KdV (mKdV) 
or Gardner (KdV-mKdV)-like equations. The properties of the phase velocity and the Landau 
damping rate of DASWs are studied for different values of the ratios of the temperatures (cr) and 
the number densities (fi) of hot and cold ions as well the cold to hot ion mass ratio m. The distinctive 
features of the decay rates of the amplitudes of the KdV, mKdV and Gardner solitons with a small 
effect of Landau damping are also studied in different parameter regimes. It is found that the 
Gardner soliton points to lower wave amplitudes than the KdV and mKdV solitons. The results 
may be useful for understanding the localization of solitary pulses and associated wave damping 
(collisionless) in laboratory and space plasmas (e.g., the F-ring of Saturn) in which the number 
density of free electrons is much smaller than that of ions and the heavy, micron seized dust grains 
are highly charged. 

PACS numbers: 52.27.Cm, 52.35.Mw, 52.35.Sb 


I. INTRODUCTION 

Dusty plasmas with an admixture of two species of ions 
are ubiquitous in both space (e.g., Earth’s D and lower E 
regions, Titan’s atmosphere, F-ring of Saturn etc. ) mm 
and laboratory plasmas EH- In addition to the main 
positive ion (light or hot), dusty plasmas also contain 
a second, micron-sized positive or negative ion (heavy 
or cold), which characterizes the appearance of an addi¬ 
tional wave mode in the plasma HU- There are in fact 
several motivations for investigating collective plasma os¬ 
cillations and nonlinear coherent structures like solitons 
and double layers in dusty plasmas with two groups of 
ions (See, e.g., Refs. SHE])- In planetary rings, in particu¬ 
lar, in the F-ring of Saturn an anomalous situation may 
arise in which massive dust grains capture almost all the 
electrons from the background plasma, i.e., the number 
density of free electrons becomes much smaller than the 
ions PE|. Such density depletion of electrons (by aerosol 
particles) have also been observed in the summer polar 
mesosphere at about 85 km altitude m as well as in the 
laboratory [T21 . On the other hand, plasmas with two 
species of positive ions have been found to have appli¬ 
cations in secondary electron emission guns (e.g., used 
in bacterial inactivation) in which helium and argon are 
used as the main gas and an auxiliary gas for ionizing 
helium atoms j!3j . 

In dusty plasmas, the size of the dust grains is, in 
general, not a constant but may vary in the range of 
0.05 — 10/irn, their mass irid ~ (10 6 — 10 12 ) rrii, where 
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rrii is the ion mass, and they have atomic numbers 
Zd ~ 10 3 — 10 5 (See,e.g., Ref. [HJ. It is well known 
that the presence of charged dust grains modifies the 
dispersion properties of plasma waves and gives rise to a 
new low-frequency wave mode, called dust-acoustic wave 
(DAW), which was first theoretically predicted by Rao 
et al. US!- These waves, associated with the motion of 
highly charged and massive dust grains, occur at rela¬ 
tively long time (~ 0.1 s) and space scales (~ 1 cm). 
Furthermore, in dusty plasmas, the presence of second 
ion species can significantly modify not only the disper¬ 
sion properties of DAWs but also some nonlinear coherent 
structures (See, e.g., Refs. EH9]). These waves can admit 
collisionless damping due to the resonance interaction of 
trapped and/or free particles with the wave, i.e., when 
the particle’s velocity becomes closer to the wave phase 
velocity [TB]. Such wave damping was first theoretically 
predicted by Landau m in 1946, and later confirmed 
experimentally by Malmberg et al. m- Motivated by 
this invention, Ott and Sudan mm first studied the 
effects of (linear) electron Landau damping on the non¬ 
linear propagation of ion-acoustic solitons in electron-ion 
plasmas. They had neglected the particle’s trapping ef¬ 
fects on the assumption that the particle trapping time 
is much longer than that of Landau damping. Their evo¬ 
lution equation was in the form of a Korteweg-de Vries 
(KdV) equation which contains, apart from the lowest 
order nonlinear terms and the term that corresponds to 
linear dispersion relation of ion-acoustic waves, a source 
term that models the lowest-order effects of resonant par¬ 
ticles. In their work, they showed that a wave steepening 
of an initial wave form may occur or not depending on 
the relative size of the nonlinearity compared to the Lan- 
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dau damping. The latter was also shown to slow down 
the wave amplitude with time. 

Recently, a number of attempts have been made to 
study the effects of Landau damping on nonlinear solitary 
waves owing to its importance in different plasma envi¬ 
ronments. To mention few, Mukherjee et al. [2D] studied 
the effects of linear Landau damping on the nonlinear 
ion-acoustic solitons in a quantum plasma with the de¬ 
scription of KdV equation. The similar investigation was 
made by Barman et al. [2T j for nonlinear dust-acoustic 
solitary waves in a dusty pair-ion plasma. In a linear 
theory, the effects of dust size distribution on the Lan¬ 
dau damping rates have also been studied by Zhang et 
al. in a multi-component dusty plasma 22]. However, 
all the above works were based on either KdV equation 
for nonlinear waves or linear dispersion laws for damp¬ 
ing rates. As is well known that the KdV solitons are 
formed for parameter values well below the critical val¬ 
ues (at which the nonlinear coefficient of the KdV equa¬ 
tion vanishes) and the modified KdV (mKdV) solitons 
are formed for parameter values close to the critical val¬ 
ues. However, when both the KdV and mKdV equations 
fail, Gardner (or KdV-mKdV) equation, which general¬ 
izes the KdV and mKdV equations, is more appropriate 
to describe the evolution of solitary waves in plasmas. 
Such (Gardner) solitons are formed for parameters very 
close to and around the critical values of the plasma pa¬ 
rameters. Though, there are a number of works in the 
literature dealing with the Gardner solitons in plasmas 
(See, e.g., Ref. [23j), however, all are based on fluid theory 
approach and thus without the effects of Landau damp¬ 
ing. 

In this work, we investigate the Landau damping 
effects on dust-acoustic (DA) solitary waves (SWs) 
(DASWs) in a dusty plasmas with two different species of 
positive ions namely heavy or cold and light or hot ions. 
Following Ott and Sudan HBJI19], we derive a KdV as well 
as mKdV and Gardner equations (with Landau damp¬ 
ing) by considering lower and higher order nonlinearities 
to describe the evolution of dust-acoustic solitary waves 
for parameter values far, near and around their critical 
values. The properties of the KdV, mKdV and Gardner 
solitons with the effects of Landau damping, the linear 
phase velocity as well as the Landau damping rate are 
also analyzed with the effects of the ratios of the number 
densities (/x) and temperatures (a) of hot to cold ions as 
well as cold to hot ion mass ratio m. It is shown that the 
landau damping has an effect of slowing down the wave 
amplitude with time, and the Gardner soliton points to 
lower wave amplitude. 


II. BASIC EQUATIONS 

We consider the nonlinear propagation of dust-acoustic 
waves (DAWs) in an unmagnetized collisionless electron 
free dusty plasma which consists of extremely massive, 
negatively charged, micron-seized cold dust fluids and 


singly charged hot (light) and cold (heavy) positive ions. 
We assume that the dust particles have constant mass 
and charge. The latter, in general, can vary and intro¬ 
duce a new low-frequency wave eigenmode as well as a 
dissipative effect (other than the wave damping) into the 
system. However, we neglect this charge fluctuation ef¬ 
fect on the assumption that the charging rate of dust 
grains is very high compared to the dust plasma oscilla¬ 
tion frequency. We also assume that the size of the dust 
grains is small compared to the average interparticle dis¬ 
tance, and the ratio of electric charge to mass of the dust 
grains remains much smaller than those of ions. At equi¬ 
librium, the overall charge neutrality condition satisfies 

ficO V n m ZdUdOi (1) 

where rij o is the unperturbed number density of species 
j (where j =c, h and d, respectively, stand for cold ions, 
hot ions and charged dust grains) and z d (> 0) is the 
unperturbed dust charge state. 

The basic equations for the dynamics of charged dust 
particles together with cold and hot ions in one space 


dimension are 

d t n d + d x (n d u d ) = 0, (2) 

d t u d + u d d x u d = -( qd/m d )d x <t> , (3) 

dtfj + vd x fj - (qj/rrij )d x <fid v fj = 0, (4) 

d x (j) = -4ire(n c + n h - z d n d ), (5) 


where d a f denotes partial differentiation of / with re¬ 
spect to a. Also, Eq. Q is for the cold (j = c) and hot 
(j = h) ions with their number densities given by 



In equations the symbols n d , u dl q d (= —z d e, e 

being the elementary charge), m d , respectively, denote 
the number density, fluid velocity, charge and mass of 
dust grains. Also, v is the particle’s velocity, f :j is the 
velocity distribution function and rrij is the mass of hot 
and cold ions. Furthermore, q.j = e, (f> is the electrostatic 
potential, and x and t are the space and time coordinates. 

In what follows, we normalize the physical quantities 
according to u d -A Ud/c s , <f> -A ecf/ksTh, rij -A rij/rijo , 
n d -A rid/rido, fj -A fjV t j/n j0 , v -A v/v th , where 
c s = \/Zdk B T h /m d = ujpdXo is the DA speed with iu pd = 
\J^ndoZ^e 2 /m d and \ D = \/ k B T h /^nn d oz d e 2 denoting, 
respectively, the dust plasma frequency and the plasma 
Debye length. Here, ks is the Boltzmann constant, Tj is 
the thermodynamic temperature and Vtj(= \fksTj / rrij ) 
is the thermal velocity of j-species particles. The space 
and time variables are normalized by L and L/c s respec¬ 
tively, where L is the characteristic scale length for varia¬ 
tions of the dependent variables rij , u d , <f>, fj etc. Thus, 
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from Eqs. we obtain the following set of equa¬ 

tions in dimensionless form: 

d t n d + d x (n d u d ) = 0, (7) 

dt^d, T U d d x U d — d x (f)^ (8) 

Sd t fj + vd x fj - (m h /m j )d x <l>d v f j = 0, (9) 

(^d/L 2 ) d 2 x (j) = -{.Hcdfic + vhdn h + n d ), (10) 

rij = ^JmjT h /m h Tj / fjdv, (11) 

J — OO 

where 5 = y/z d mh/m d and p, j d = rijo/z d n d Q is the den¬ 
sity ratio for j = c, h. The charge neutrality condition 
([I]) reduces to 


and for L Ad (valid for long-wavelength perturba¬ 
tions), Eqs. ([7])-( 10) yield the following linear dispersion 
law 12T1 


v p = u}/k = (n hd + an c d) 1/2 , 


(16) 


where ui (k) is the wave frequency (number) of the plane 
wave perturbations and er = is the ratio of hot to 

cold ion temperatures. Equation ( |16| ) has the same form 
as Eq. (15) in Ref. I^Tl obtained in dusty plasmas with 
positive and negative ions. From Eq. (16) it is seen that 


the DAWs become dispersionless in the long-wavelength 
oscillations with the constant phase speed v p (< c s ). This 
implies that the time derivatives of all physical quantities 
should vanish in a frame moving at the speed v p , and one 
can thus expect slow variations of the wave amplitude in 
it. So, we introduce the stretched coordinates as [22 


£ = e 1/2 (x — Mt), t = e 3/2 t, 


(17) 


dcd + dhd = 1, or Hhd = n /(1 + n), 


( 12 ) 


where fi = riho/n c o is the ratio of hot to cold unperturbed 
ion number densities. As in Ref. [22 there are three basic 
parameters. The first one is the ratio 6 = \Jz d mh/m d 
or md (m = m c /mh), which represents the finite inertial 
effects of hot and cold ions, and, in particular, the Lan¬ 
dau damping. The second one is the ratio of perturbed 
dust density to its equilibrium value n d i/n d o, which mea¬ 
sures the strength of nonlinearity in electrostatic per¬ 
turbations. Lastly, the parameter Ad/A, measuring the 
strength of the wave dispersion due to deviation from 
the quasineutrality. This parameter will eventually dis¬ 
appear in the left-hand side of Eq. (101, if one considers 


the normalization of £ by Ad instead of L. Note that 
since the dust particles are assumed to be cold, the Lan¬ 
dau damping effects are provided solely by the hot and 
cold ions, and the damping rate is oc 5. 


III. DERIVATION OF KDV EQUATION 


Before going to derive the KdV equation we first con¬ 
sider the following scaling for the basic parameters de¬ 
scribed in the previous section as im nu 

d = a ie, n dl /n d0 = a 2 e, A 2 D /L 2 =a 3 e, (13) 


where e (< 1) is a smallness positive parameter and oq 


o(l) (j = 1,2,3) is a constant. Under the scaling (13) 
Eq. § for hot and cold ions reduces to 


ceiedtfh + vd x f h - d x <f>d v f h = 0, 


(14) 


where M is the nonlinear wave speed (relative to the 
frame with x and t as coordinates) normalized by c s , to 
be shown to be equal to v p later. The dynamical variables 
are expanded as 


n d = 1 + aztri'P + a 2 e 2 n^ 2) -|-, 

u d = a 2 eu { -p + ale 2 u { p H-, 

4>= a 2 e<j>^ + a^e 2 d-, 



1 + a 2 enf' 1 + a 2 e 2 nj 2) -|- 

/j 0) + a 2 e/j (1) + oi\ e 2 /j 2) + 


(18) 


where the equilibrium distribution function /j°\ j = 
c, h , are assumed to be the Maxwellian 


/j 0) = \/lj2 rrexp [(mjT e /m e Tj)(-v 2 /2)\ . (19) 


In the next step we substitute the expressions from 
Eqs. (17) and (18) into Eqs. Q, 0, @ 

and (15), and equate different powers of e. In this way, 
we obtain the equations for first and second order per¬ 
turbations. Since the method and intermediate steps for 
the derivation of the KdV equation are almost similar 
to that in Ref. H2 we skip presenting first and second 
order expressions for the perturbations. However, those 
can be given for the derivation of the Gardner equation. 
The mKdV equation will then directly follow from the 
Gardner equation. Thus, following Ref. [22 we obtain 
the following KdV equation 


dn „ [°° dn d£' 


, dn d 3 n 

bn iK +c W = ’ 


( 20 ) 


and 


oiiedtfc + vd x f c - (1 /m)d x (j>d v f c = 0. (15) 

Next, in the small-amplitude limit (i.e., e —> 0) in 
which inertial effects of hot and cold ions are neglected, 


where the coefficients of the Landau damping (a), the 
nonlinear ( b ) and the dispersive term (c), respectively, 
are given by 




(y/m-cT 3/2 )nhd 


( 21 ) 
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b = (a 2 /2yfaa{) (3 - <j 2 /<jI) . ( 22 ) 

c= (a 3 /2)(aa 1 )~ 3/2 , (23) 

with <ti = 1 + (cr -1 - 1 )n h d = At c d(cr + n)/a, a 2 = 
1 + (ct ~ 2 - 1 )/jm = tt c d(o- 2 + aO/ct 2 and ct 3 = 1 + (a~ 3 - 
l)^i^d = fi c d(cr 3 + n)/<J 3 - We note that the nonlin¬ 
ear coefficient b of the KdV equation ( f20j ) vanishes for 
Ri = 3- a 2 /a\ = 3- (1 + ^){a 2 + /j,) / (a +/j) 2 = 0 along 
a curve in the a — /r plane (See Fig. [TJ . In this case, the 
KdV equation fails to describe the nonlinear evolution 
of DASWs, and one has to look for some higher order 
perturbations, i.e., the evolution of DASWs can then be 
described by rnKdV or Gardner equations to be obtained 
shortly. Furthermore, in absence of the Landau damping 



FIG. 1. Contour plots of Ri = 3 — a 2 /cr 2 = 3 — (1 + 
/i)(cr 2 + + fj,) 2 = 0 (The thin or black curve) [i.e., van¬ 

ishing of the nonlinear coefficients b and bi of the KdV and 
Gardner equations (20 1 and (871] and R 2 = 1 — 03 /15<ti = 
l — (l+IJ,) 2 (<j 3 +[i)/(a+IJ ,) 3 = 0 (The thick or red curve) [i.e., 
vanishing of the nonlinear coefficient b 2 of the mKdV equation 
( |90| )] showing the critical regions in the a — /r plane. There 
is a common critical point (cr, /r) « ( 10 , 10 ) at which both the 
KdV and mKdV equations fail to describe the evolution of 
DAWs. 


the KdV equation admits compressive or rarefactive soli- 
ton solutions according to when b > 0 (b < 0 ), i.e., when 
b falls in the region i?i > 0 (i?i < 0). We will, however, 
later see that the region for which b < 0 is inadmissible 
for the existence of KdV solitons with a small effect of 


the Landau damping [See Eq. (1201 in Sec. VI]. We also 
mention that the KdV solitons with Landau damping are 
formed for values of the parameters cr and /1 far below the 
critical values along the curve R\ = 0, the mKdv solitons 
are formed for the parameters close to the critical values 
and Gardner solitons are formed approximately at and 
around the critical values of cr and /i. Thus, the KdV 
equation may be valid for the parameters <r and /r satis¬ 
fying Ri ~ o(l). However, in the case of R\ ~ o(e) or 
Ri ~ 0, one has to deal with mKdV or Gardner equations 
which involve higher-order (of e) perturbations for the 
evolution of dust-acoustic solitary waves in the plasma. 

Figure [l] shows the contour plot of Ri = 0 (the thin or 
black curve) in the a — fi plane in which the nonlinear co¬ 


efficients b and b \, respectively, of the KdV and Gardner 
(obtained later) equations (20 1 and (87) vanish. The plot 
of R 2 = 0 (The thick or red curve) represents the same, 
however, for the vanishing of the nonlinear coefficient b 2 
of the mKdV equation (90) (shown later). We find that 
when cr and p, assume values satisfying R\ ss 0 the KdV 
solitons may not exist, the evolution of dust-acoustic soli¬ 
tons can then be described by the mKdV equation. The 
latter can also be inappropriate for the evolution of soli¬ 
tons when R 2 ~ 0 for some values of a and /i. From Fig. 
[Tj it is also seen that the curves for Ri = 0 and R 2 = 0 
has a common point of intersection (cr, /r) « ( 10 , 10 ) at 
which both the KdV and rnKdV solitons cease to exist. 
Also, the lower parts of these curves below the intersec¬ 
tion point are close to each other, implying that the pa¬ 
rameters close to the lower part of R\ = 0 may be some 
what close to that of R 2 = 0. In this case, the Gardner 
equation is necessary for the evolution of DASWs, which 
will be shown later. 


IV. DERIVATION OF GARDNER AND MKDV 
EQUATIONS 


To derive the Gardner equation for DASWs we con¬ 
sider Ri ~ o(e). In this case, the scaling for the basic 
parameters [Eq. (13)] and the equations for hot and cold 
ions, i.e., Eqs. (141 and (15) are to be changed by the 
substitution e ->r as we are looking for perturbations 
with higher-order (of e) effects. Thus, we have the mod¬ 
ified scaling for the basic parameters 

5 = aqe 2 , n dl /n d0 = a 2 e 2 , A 2 D /L 2 = a 3 e 2 . (24) 


The modified kinetic equations for ions are 

otie 2 d t f h + vd x f h - d x (j>d v f h = 0, (25) 

and 

aie 2 d t f c + vd x f c - (1 /m)d x (j}d v f c = 0 . (26) 

Also, the stretched coordinates are redefined as 

£ = e(a: — Mt), t = e 3 t, (27) 


However, the dynamical variables are expanded in the 
same way as Eq. (18) for the KdV equation. Next, we 
substitute the expressions from Eqs. (18) and (27) into 
Eqs. 0 , ©, jnj, ( f25| ) and ( [26] ) and equate differ¬ 

ent powers of e. We temporarily drop the constant a 2 in 
the subsequent expressions and equations. However, we 
note that the constants aq, a 2 and a 3 will explicitly ap¬ 
pear in the coefficients of the Landau damping, nonlinear 
and the dispersive terms in the evolution equation. The 
results will be given in the following subsections. 


A. First-order perturbations and nonlinear wave 
speed 


We successively equate the coefficients of e 2 from Eqs. 
0 and 0, the coefficients of e from Eqs. (10) and ( |Tl| ), 
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and the coefficients of e 2 from Eqs. (25) and (26) to 
obtain 


= «dV M . 

0 = -(fleeing + Hhdn^) + n^\ 

,- /»oo 

= ^JmjT h /m h Tj / ff^dv, 

J — OO 

^/h 1} + v fh )d ^ {1) = °> 

+ erv /i 0) «^ (1) = o. 

Equations (p8| and (p9|) reduce to 


n ? 


Ai) 


= -<j>^/M 2 . 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


Also, Eq. (32) yields the following relation 


dsfh } = ~fh ) d i (l) {1) + A(£, t)S(v), (35) 

where S(v) is the Dirac delta function and A(£, r) is an ar¬ 
bitrary function of £ and t implying that the solution for 
df}P/d£ is not unique. Thus, for the unique solution to 
exist we add a higher-order term e 6 ai (df^/dr^ origi¬ 
nating from the term e 5 ai ( dfh/dr) in Eq. ( |25| after the 
expressions (27) and ((l8| being substituted. Thus, we 
rewrite Eq. (32) as 


<?d T fhJ + vd^ = -vf^ ] d^ w , (36) 

Similarly, from Eq. ( |33[ ), we obtain 

one 2 d T f^ + vdzf£ ] = -ctu/ c (0) <9 5 </ (1) . (37) 

and (f37l) can be given 


The unique solutions of Eqs. 
as 


fj 1] = Um /j e 1} , 

J e->0 J 


(38) 


where (j = h, c) are to be obtained by taking the 


Fourier transform of Eq. (36) with respect to £ and r 
according to the formula 


/ oo /*oo 

/ /(£,r)e^—)d£rf r . 

-OO -J — oo 

Thus, we obtain 

fhJ = - kv fh ] ( kv ~ e 2 «iw) 1 0 (1) . 


(39) 


(40) 


Next, to remove the singularity in Eq. (40j) we replace w 
by u> + ir] with r] (> 0) being small, and obtain 


Proceeding to the limit as e —► 0, using the Plemelj’s 
formula (in which P denotes the Cauchy principal value) 

lim (x + ie)~ 1 = —iTrS(x) + P (1/a;), (42) 

e—^0 

and noting that xP(l/x) = 1, xS(x) = 0, we obtain 

= -/W 1} - (43) 


Next, we take the Fourier inversion of Eq. (43) to obtain 

(44) 


if = -if f 


>h ~ Jh 

Similar expression for cold ions can also be obtained from 


Eq. (37) as 


if = 


(45) 


Substituting the expressions from Eqs. (44) and (45) into 
Eq. (31), we obtain 


4 1} = -f 


(46) 


n f = — of 1 ). (47) 

The expressions for nf (j = c, h, d) from Eqs. (34), (46) 


and (47) can now be used in Eq. (30), to obtain the 


following expression for the nonlinear wave speed [21] 

M = (n hd , + <JHcd)~ 1/2 = Vi 1 + f/f + o-). (48) 

Evidently, M has exactly the same expression as v p ob¬ 
tained in the linear dispersion law ( |16| ) for plane wave 
perturbations. Figure [2] shows the characteristics of M 
with a for different values of the density ratio fi. We 
note that since a = T//T c > 1, the quantity under the 
square root in Eq. (48) is always smaller than unity, i.e., 


the DAWs propagate with a phase speed being smaller 
than the DA speed c s . We also find that as the value of 
H increases, the value of the phase velocity M increases. 
However, the same decreases for increasing values of the 
temperature ratio a, and approaches to zero for a ^$> 1. 
The latter is, however, inadmissible in the present theory 
as it would not give any physical result. We mention that 
the results in this subsection (up to the linear theory) are 
the same as for the KdV and Gardner or mKdV equa¬ 
tions. However, the results (nonlinear) will be different 
in the expressions for second and higher order perturba¬ 
tions. 


B. Second-order perturbations 


We successively equate the coefficients of e 3 from Eqs. 
([7]) and j8]), the coefficients of e 2 from Eqs. (10) and (11) 


as well as the coefficients of e 3 from Eqs. (141 and (15) 
to obtain 


/f = [(kv - e 2 aiui) - irjaie 2 ] 1 <f(41) 


f (0) 


-Mnf + uf = -nf uf, 


(49) 
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FIG. 2. The linear phase velocity M [Eq. (481] of DAWs is 
shown with respect to the temperature ratio a (= Th/T c ) for 
different values of the density ratio fi (= nho/rico) as in the 
figure. 


u d ] = - (1 /M)^ + (1/2M 3 ) (V (1) Y 


(50) 


( 2 ) ( 2 ) , ( 2 ) 
^d — l^cd'^'c M hdTlfo ? 


(51) 


( 2 ) 

n 3 = 


J - poo 

\JmjT h /m h Tj / /j 2) cfo, (52) 

J — OO 


vd dh 2) - 9^ l) dvfh ] + = 0, (53) 

vdzfP - (1 MdtfWdyfP + avfWdtfW = 0. (54) 
Eliminating u^ from Eqs. (49) and (|50|) and using the 


first-order expressions from Eqs. (29) and (34), we obtain 

2 


/ 2 ) 

l d 


(3/2M 4 ) (<A (1) )" - (1/M 2 )<Y 2) . (55) 


Substituting the expressions for fj from Eqs. (44) and 
(45) into Eqs. ([53]) and ([54]) we successively obtain 


vd ifh ] ~ ^/YV (1) 9{</> (1) + v d^ 2) = 0, (56) 


vd^f^ — cr 2 u/ c (0) </> (1) ^/> (1) + crvf^dtcj)^ = 0 . (57) 

As in the previous section, to obtain the unique 
( 2 ) 

solutions for /> ' we introduce a higher-order term 
e 7 ai {df^ /dr^j originating from the term e 5 ai ( dfj/dr) 
in Eqs. [25| and ([26]) after the expressions ([27| and (fT 8 [) 

^ ( 2 ) (2) * 1—1 

being substituted, and replace /> ' by /> '. Thus, we 


rewrite Eqs. (56) and (|57|) as 


“1 e2 drfL +vd i fi ) +vfl ] d ^ 
= (l/2)t,/<°>0 c (0«) 2 , 


( 58 ) 


ai e 2 d T fW + vdsfW + avfWd&W 
= (a 2 /2 )vf^dd<t> W ) 2 - 


(59) 


( 2 ) 

As before, the solutions /> , j = h,c can be found 
uniquely from Eqs. (58) and (59) by the following re¬ 
lation 


f? ] = lin4 /i 2) ‘ 

J e->0 J 


(60) 


So, taking the Fourier transform of Eq. (58) according 
to the formula (391, and replacing oj by u + ir) (77 > 0 ), 
we obtain 


A 2) _ 

J he 


-a>+)»“')= 


k «fh°’ 


(kv — e 2 ai w) — ir]a\e 2 

( 61 ) 

Proceeding to the limit as e —> 0 and using the Plemelj’s 
formula (42) together with xP(l/x) = 1 and xS(x) = 0 , 
we obtain from Eq. © the following 


A 2 ) _ 

Jh — 


-Y 2 ) + (l/2 ){^f 


AO) 
Jh ■ 


(62) 


Next, the Fourier inverse transform of Eq. (62) yields the 
following solution for the hot ion distribution function 

-Y 2 ) + (l/2)(</W) 2 l fl 0) . (63) 


A2) _ 
Jh ~ 


By the same way as above, one can also obtain the similar 
expression for the cold ion distribution as 


/i 2> = 


-offP^ + (a 2 / 2) (V 1 ' 1 ) 


/j 0) . (64) 


The expressions for the second order perturbations of the 
number densities can now be obtained from Eq. (52) by 
using Eqs. (63) and ( |64| ) as 

n h } = -</> (2) + (1/2) (</> (1) ) 2 , 


and 


+ (a 2 /2)(Y 1} ) 2 


(65) 


( 66 ) 

Furthermore, substituting the expressions for rS 2 \ n^\ 

n^ from Eqs. (55), (65) and ( 66 ) into Eq. (51), we 
obtain for nonzero perturbations the following relation 

— (1/2) (vhd + AGdO- 2 ) (/> (1) ) 2 + (Hhd + Acdv) (/ (2) 

= (1 /Af 2 )</> (2) - (3/2M 4 )(</ (1) ) 2 , (67) 

in which the coefficient of <// 2 ) vanishes by the dispersion 
relation (48) yielding 


b = (l/2)(crcri) 2 (3 - 0 - 2 / 0 - 2 ) 


m = 


= 0 , 


( 68 ) 


where cy = (cr+/j,)/cr(l+n) and a 2 = (cr 2 + /z)/cx 2 (l + /z). 
Here, we note that since is of the order of e, the 
factor i?i = 3 — er 2 /cr 2 must be at least of the order of 
e, which we have already assumed in the beginning of 
this section. So, the term b will contribute to the terms 
in the coefficient of e 3 of the Poisson’s equation in Sec. 
IVC [See Eq. ©], and finally in the coefficient b\ of the 
Gardner equation (f87[. 
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C. Third-order perturbations 


Equating the coefficients of e 4 from Eqs. (pf] ) and ([8]), 
the coefficients of e 3 from Eqs. (10) and (11), and the 
coefficients of e 4 from Eqs. (25) and (26), we successively 
obtain 


-Mdi-nf + d T n { d 1] + d^P + d^npuP) 
+ d i( n d )u d ) ) = °> ( 69 ) 


druP - M<%4 3) + di(u { d ] u { d ] ) = <%</> (3) . (70) 


dl4> {1) =-(nhdnP + He dnP) + n ( d ] 

+ (l/2)(crcri) 2 (3 - cr 2 /crj) (0 (1) ) , (71) 


I - poo 

Pp = nrijT h /m h Tj / ff^dv. (72) 

J — OO 


-aiMd^fP + vd^fP + vfPd^ 3) 

-d^ 2) d v fP - d^d v fl 2) = 0, (73) 


-aiMd^fP + vd^fP + avfPd^p 3 ^ 
-(1/m) ( d^d v fP + d^d v fP) = 0. 


(o\ 

As before, to get the unique solutions for /' (for j = 
h , c), we introduce a higher-order term e 8 «i (dfP /d 
originating from the term e 5 ai (dfj/dr) in Eqs. 

(26), and then replacing fp by fP, where 


(74) 


T 

and 


be 


f (3) = lim f^ 3) 


we obtain from Eqs. (73) and (74) as 

aie 2 d T fP +vdtfP +vfPdi:p 3) 
= (D ha + vD hb )fp , 


(75) 


(76) 


aie 2 d T fP + vdsfP + avfPdtfW 
= (D ca + vD cb )fP, (77) 


where we have used Eqs. ([44]) and (63) to obtain Eq. 
(76), and Eqs. (451 and (|64|) to obtain Eq. ( [77] ). The 
expressions for D ba , Dhb , D ca and D cb are given by 


D ha = 

D hb = -(l/6)<%(</ (1) ) 3 + dc (^ 4 V (2) ) , (78) 

D ca = -a 1 Mad^p 1) , 

D c b = — (cr 3 /6)c? € (</> (1) ) 3 + er 2 <% ■ 


Next, taking the Fourier transforms of Eq. (76) with 
respect to £ and r according to the formula (39) and 
replacing ui by to + ir) [rj > 0), we obtain 


7(3) _ _ 
J he 


kvfpP 3 ') + i{D ha + vD hb )fp 


(kv — e 2 aiw) — ir/aie 2 


(79) 


In the limit of e —> 0, Eq. (79) yields 


fh + fh = — * P iX/kv) + J7r(sgn k/k)S(v)] x 


(D ha + vD hb )fP, 


(80) 


where we have used the Plemelj ’s formula (42) together 
with ccP(l/a;) = 1, xS(x) = 0 and S(kv) = (sgn k/k ) S(y). 
Now, multiplying both sides of Eq. ( |80| by ik and then 
integrating with respect to v, we obtain 


ik (nP + 


= D hb + sgn (k)D ha 


(81) 


which under the Fourier inverse transform reduces to 

d(_nP +d£<p^ =d((<j>^(f)^) - (l/6)c?^(</> ( ^ 1 ^) 3 

+\frp 2 F _1 (i sgn (k)D h Pj . (82) 

Here, F^ 1 stands for the Fourier inverse transform. Us¬ 
ing the convolution theorem of Fourier transform and 
noting that F -1 [i sgn(fc)] = — (l/-7r)P (l/£), we obtain 
from Eq. ( |82| ) 

dt n P +<%</> (3) = 9 5 (</> (1 V (2) ) - (l/6)9j(^ (1) ) 3 

+cnM (l/v^J P / %/> (1) ^ 77 - (83) 

Similar equation for the third-order perturbations can be 
obtained for cold ions as 


(l/i/mcr) d^nP + 





d^(jX x ^ 





(84) 


D. Gardner and mKdV equations with Landau 
damping 


In order to obtain the evolution (Gardner) equation 
for first-order perturbations, we first eliminate uP and 
np from Eqs. to obtain 


dfp 1] = -Hhdd^nP - Hcdd^nP - (1 /M 2 )d 5 </> (3) 

+ (1/M) d T nP + (1/M) (d T uP + <% 


"d “d 


( 1 ) ( 2 ) ( 2 ) ( 1 ) 
n ^uy + ny u y 


+P ( f 

+ ((J(Ti) 2 (3 - cr 2 /cr 2 ) ^d^(p V> . 


( 85 ) 


































Next, we use Eqs. (831 and (84) to eliminate 
n^ from Eq. (|S5|) , and obtain 


(3) 


and 


<9f<(> (1) = - (fi hd + a 2 fi cd ) <% 
+( 1 / 6 ) (n hd + cr 3 fi cd ) (/> (1) ) 


a\M 


f-^hd 


mcr 3/2 n c d 


) p r ^ 

' J — oo 


drU^ 


n d Ur, +n d u d 


(1/M) (c 




' —OO 

,( 1 ) 




n ( (1) (2) 

( u\ V 


*d u d 


( 86 ) 


V2n V 

+(1 /M) 

+d e (r (1) ~ (2) 

+ (trcr 1 ) 2 (3 - o- 2 /a 2 ) 

( 2 ) ( 2 ) 

Next, we substitute the expressions of u d and n d from 
Eqs. (50) and ( |55| into Eq. ( 86 ). Here, we note that 
the coefficient of is ~ (3 — a 2 /crf) ~ e, and so 

the term oc (3 — < 7 2 /cr/) (<// 2 ) <j)^) ~ e 4 [ since 0 ( 2 ) ~ e 2 
and < fit 1 ' 1 ~ e ]. Thus, this term will not contribute to 
Eq. ( 86 ) and we discard it. Finally, we replace by 

and 


TO 


using Eqs. 


to obtain from Eq. 


( 86 ) the following Gardner (KdV-mKdV) equation 
dn „ dn d£ 

a? + 


, dn 2 dn d 3 n 

+b ' n si + lKn ai +c W = 


(87) 


where the coefficients a and c, corresponding to the Lan¬ 
dau damping and dispersive effects, are the same as for 
the KdV equation ( |20| ) . The nonlinear coefficients b\ and 
62 are, however, different and given by 


h = (3 - a 2 /af) 


( 88 ) 


(15/4) (a 2 /(crcri) 1/2 ) [1 - (<7 3 /15erf)] . (89) 

Here, the expressions for <j\ and cr 2 are given before and 
cr 3 = (a 3 + fi)/a 3 (l + fi). If we neglect the strength of 
the Landau damping i.e., if we set oq = 0, the term oc a 
in Eq. (87) vanishes, and Eq. (87) reduces to the usual 


b , = 


KdV-mKdV or Gardner equation for the propagation of 
DASWs in dusty bi-ion plasmas. Equation (87) is the 


required evolution equation for the first-order dust den¬ 
sity perturbations when the plasma parameters <7 and /i 
assume values close to the critical values, i.e., Ri ~ o(e). 
In fact, this equation is well applicable to DASWs when 
the plasma parameters assume the critical values or are 
very close to or around the critical values of a and fj,. 

On the other hand, when i?i = 0, i.e., when the param¬ 
eters cr and n assume the critical values along the thin (or 
black) curve as in Fig. [lj the evolution of DASWs will 
no longer be describable by the KdV equation ( p0| , but 
the mKdV equation, which can be obtained by setting 
6i=0 in Eq. (871 as 


dn 

dr 


aP 


f 00 dn dg 


b 2 n 


,dn 

di 


d 3 n 

c w =0, 


(90) 


where the coefficients a, b 2 and c are the same as in 
the Gardner equation (87). We remind here that the 


evolution of DASWs with the effects of Landau damp¬ 
ing are governed by the KdV, rnKdV or Gardner equa¬ 
tions depending on whether the density ratio /r and the 
temperature ratio cr assume their values well below the 
critical values, very close to the critical values or near, 
around and exactly the critical values of /i and a along 
the curve in Fig. |T] Note that the mKdV equation is ob¬ 
tained on the assumption that R± = 0. However, when 
R 2 = 1 — < 73 / 1517 / = 0, i.e., the nonlinear term b 2 in the 
mKdV equation vanishes in the <7 — fi plane (See the thick 
or red curve in Fig. [l]), then the mKdV equation fails to 
describe the evolution of DASWs. From Fig. [l] it is also 
clear that there are some values of cr and /i which may 
be closer to both the curves (e.g., the lower parts of the 
curves for i?i = 0 and R 2 = 0 below the point of inter¬ 
section). In these parameter regimes both the KdV and 
mKdV equations may not be appropriate for the descrip¬ 
tion of DASWs in plasmas. Thus, the Gardner equation, 
being the generalized form of KdV and mKdV, describes 
the evolution of DASWs in dusty bi-ion plasmas for a 
wide range of values of <7 and fi including the critical 
values as in Fig. [l] 


V. LANDAU DAMPING RATE 

To obtain the Landau damping rate for DASWs we 
take the Fourier transform of Eq. (871 with bi , b 2 and 
c being set to zero, according to the formula (39), and 


use the result that the inverse transform of i sgn(fc) is 
— (1/ 7 r) P(l/£). Thus, we obtain the Landau damping 
rate as |2L |y| = na. This damping rate is shown in 
Fig. 1 with respect to the temperature ratio cr, and for 
different values of the density ratio fi and the mass ratio 
m. We find that for dusty plasmas with equal mass of 
positive ions, the value of |y| decreases (increases) with 
an increasing value of the parameter <7 (fi). However, the 
case with unequal mass of ions (m > 1 ) is quite distinc¬ 
tive. In this case, the damping rate initially increases 
as long as 1 < cr < 5, and for cr > 5, the behaviors 
of I 7 I remain similar to the case of equal mass of ions. 
Thus, in contrast to dusty pair-ion plasmas [2Tj or dusty 
electron-ion plasmas, an enhancement of the hot ion con¬ 
centration (or a reduction of the number of cold ions into 
the dust grain surface in order to maintain the charge 
neutrality), increases the Landau damping rate in dusty 
bi-ion plasmas. 


VI. SOLITONS WITH LANDAU DAMPING 

We first obtain an approximate analytic soliton so¬ 
lution of the Gardner equation (87) with a small ef¬ 
fect of the Landau damping. Here, we assume that 
1 ~ bi, b 2 ~ c a e. This can easily be verified by 


considering dusty plasma parameters satisfying < 7 , fi > 1. 
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FIG. 3. The linear Landau damping rate |-y| = na is plotted 
against the temperature ratio <r for different values of the 
mass and the density ratios m and as in the figure. 


Now, in absence of the Landau damping Eq. (871 reduces 
to 


d T >n + b/nd^n + b' 2 n 2 d^n + 9|n = 0, 


(91) 


where we have rescaled the time variable r and the co¬ 
efficients as r' = re, b\ = b\/c, b' 2 = b 2 /c. A traveling 
wave solution of Eq. (911 can be obtained as PS] 




n (£, T ') = 6C7 b[ ± Acosh ((£ — Ut')/W 


(92) 


where W = y6b' 2 /(N — b' 2 ) is the constant width and 
— —2 

U = (A — bf)/6b 2 is the constant phase speed (normal¬ 
ized by c s ) of the DASWs. It can easily be verified that 
Eq. ( |87| ) conserves the total number of particles and an 
initial perturbation of the form (92) decays to zero (See 
for some details, e.g., Ref. ED. The latter implies that 
the wave amplitude A is not a constant but decreases 
slowly with time. Thus, we introduce a new space co¬ 
ordinate in a frame moving with the solitary wave and 
normalized to its width as PH 


In what follows, we investigate the solution of Eq. (941, 


and introduce a multiple time scale analysis with respect 
to a! . Thus, we write PUPS] 

n(z, r') = n (0) + a'n {1) + a' 2 n {2) + a'V 3) + • • • , (95) 

where n^\ i = 0,1,2,3, ••• are functions of 

r' = to, ri, T 2 , 73 , • • • in which t* = a' l r'. So, 
we have dn^/dr’ = [dn^/dr 0 ){dr 0 /dr') + 
{dn^/dr 1 ){dr 1 /dr') + dn^’/dr' 

(dn^/dr 0 )(dr 0 /dr') + {dn^/dr^drjdr') + 

etc. and the similar expression for dN/dr. Substitut¬ 
ing the expansion (95) into Eq. ( f94j ) and using the 
expressions for dn^/dr' , dn^/dr' and dN/dr' we 
obtain 

(d T >n + a'd Tl n ^ H-^ + a' (d T /n^ + a'd Tl n^ H— • ^ 

+ [—A + Bz ( d T ’N + a'd Tl N +•••)] x 

d z (n (0) + «'n (1) + •••) + (l /Wj (b[ + &' 2 n (0) ) n (0) <9 2 rc (0) 

+= [b[d z (nWn(°))+^( 0 > (r»(°)a z n« + 2n< 1 >0 z n<°>) 

+ ( 1/W^ d 3 n {0) + a' (l/IT 3 ) d 3 n {1) 


1 

+a' = P 

W J -oo 




dz' 


= o, 


z — z' 


(96) 


Equating the coefficients of the zeroth and first-order of 
a', we successively obtain from Eq. (96) as 


P 


d T > + (Bz)d T 'Nd z 


n (0) + M 0 <9 2 n (0) = 0 


(97) 


P 


dr' + (Bz)dr'Nd z ] n (1) + d z M 0 n {1) = j3Rn {0 \ (98) 


—3 

where ft = W and 


M 0 = d 2 + 


i (0) (b[ + b' 2 n(°A (AW 


- 1 


(99) 


* = (f+ w/ - l N ‘ dT )' 

(93) 

where N = N(a, r') is assumed to vary slowly with time 
and n = n{z, r'). Under this transformation Eq. (87) 
reduces to 


d T m - 


—A + Bz (dN/dr' 


d z n + 


(b'JW) 


nd z n 


b' 2 /W^ n 2 d z n + (l/IR 3 ) 
(a 1 /Wj P 


d 3 r 


(j y! 

din -- = 0, 


z — z 


(94) 


where A = (A — b' 2 )/6b 2 W, B =N/ (A 2 — b' 2 ), a' = a/c 
and we have used dn/dz' = dn/dz at 2 = z'. 


i?n (0) = - |d Tl n (0) + ( Bz)d Tl Nd z n (0) 


-=P 
W J -00 


cGz (0) 


dz' 


z — z' 


( 100 ) 


Using the boundary conditions, namely, n^°\ dn^/dz, 
d 2 n^/dz 2 , d 3 n ( '°' > /dz 3 —» 0 as z —> ±00 it can be 
shown that n^ 0 -* = 6/7o/(6 , 1 ± Acoshz) is a soliton so¬ 
lution of the equation Mq (dn^/dz) = 0. Thus, = 


6Uo/(6 , 1 ± Acoshz) will be a soliton solution of Eq. (97) 
if and only if dN/dr = 0. So, Eq. ([98]) reduces to 


P d T n^ + d z M 0 n^ = j3Rn^°K ( 101 ) 


Next, for the existence of a solution of Eq. (101), Rn ^ 
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(a) (b) 



(c) 


FIG. 4. The solitary wave amplitudes N(a,r) (normalized to their equilibrium values No) with the effect of Landau damping 
are shown for KdV, mKdV and Gardner solitons [Eqs. (1191, (|122|>, (1 121]. Subfigure 4(a) shows that when g and a assume 


values (m = a = 3, g = 7.4) well below the critical values on the curve Ri = 0 and fall in the region Ri > 0, the decay rates 
of the amplitudes of the KdV and Gardner solitons with respect to r are approximately the same as expected. However, the 
rate for the mKdV soliton is quite higher than the KdV and Gardner solitons. For values of g and a (m = 3, cr = 15, g = 4.5) 
close to the critical values, but in the regions Ri, Ro < 0 for which b, b± < 0 (and hence KdV equation fails to govern the 
dynamics of dust-acoustic solitons), subfi gure |4(b) shows that the gap between the decay rates of the Gardner and mKdV 
solitons becomes lower than that in Fig. |4(a)| On the other hand, when g or a values (m = 3, <r = 11, g = 5.3) are very 
close to the critical values, but still in the regions Ri, R 2 > 0 for which b, 61 « 0 and 62 > 0, subfigure 4(c) shows that the 
gaps among the amplitudes of the KdV, mKdV and Gardner solitons again increase, and Gardner equation seems to be more 
appropriate for the description of dust-acoustic solitons. 


must be orthogonal to all solutions g(z) of L + [g\ = 0 
which satisfy c/(±oo) = 0, where L + is the operator ad¬ 
joint to L = ( d/dz)M , defined by 


by 


L+ 



6N(N ± b\ cosh z) 
(b[ ± cosh 2 ) 2 



(103) 


ipi(z)L[i/j 2 (z)\dz = / i/j 2 {z)L + [ip 1 (z)]dz, (102) 


where iJji(±oo) = ^2 (± 00 ) = 0, the operator L + is given 


and the only solution of L + [g\ = 0 with <?(±oo) = 0 is 
g{z) = (N — b'i)/N(b[ ±7Vcoshz). Thus, we obtain 


/ °° 1 —2 — 

= (N — 6 , 1 2 )(6 / 1 ± N cosh z)- 1 Rn^ dz = 0. (104) 

-00 N _ 

Subs titut ing the expression for from Eq. (100) into 

Eq. (1041, we get 
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i- 


dN r 

' dn L 


(TV — b' 2 ) TVcosh ^ 


Nz sinh z 


N 2 ( b / iTVcosh z) 3 (6/ iTVcosh z) 3 


dz 


NW 


r°° r°° 1 d 

1 

dzdz' 

/- 00 J -00 (6/ ± TV cosh z) dz' 

(6/ ± TV cosh z') 

z — z' 


= 0, 


(105) 


r 


from which we obtain after few steps the following 







dN 

(3TV — 2b' 2 ) 

8 ^ 

to 

1 

-ts 

8 ^ 

dz 

dr 4 

2TV 2 

7-oo (6/ iTVcosh z) 2 

_/y 2 7-oo (6/ iTVcosh z) 3 

(W 

1 

'•OO nOO ^ Q 

1 

dzdz' n 

Ny/w; ^ J 

-00 7-oo (6/ ± TV cosh z) dz' 

(6/ iTVcosh z') 

z — z' 


dz 


(106) 


r 


(b[ i TVcosh *) 2 dz = 


2/(IV 2 - b' 2 ) 3/2 


Now, we perform the integrations as follows 


— , — 2 


26/tan -1 (6/ T TV)/-\/a - 6/ 2 + 


sJn 2 - b' 2 


(107) 


(b[ zbiVcosh z) 3 dz = 


2/(TV 2 - bf) 5 / 2 


—2 n 


(2&i + iV )tan 


-l 


(6'iTiV)/ViV 2 -6'i 2 ) i^ 


(108) 


Finally, using th e re sults in Eqs. (1071 and (108), we 
obtain from Eq. (106) as 


<9TV , „ f 00 1 

w—Ai — A2P / / -=- 

J —OO J — OO (6/ i TV cosh z) 


d' z (b[ i AT cosh 2 '^ [cfe'/li — /)] cfe = 0, (109) 

where Ai and A 2 are given by 

Ai = [l/(TV 2 - 6/ 2 ) 3/2 l 86/tan -1 f (6/ ^N)/y/p? - 6/ 2 


+ ((3TV + 6/ 2 )/Tv) ^ - 6/ 2 


Ao — 


(TV 2 - b' 2 ) 3/2 /Ny/W 2 


( 110 ) 


( 111 ) 


Equation (109) is a first-order (in r) partial differential 
equation for the solitary wave amplitude TV(a, r) whose 
exact analytic solution is much complicated. However, 
one can consider some particular cases as follows: 

We note that for small-amplitude perturbations 
e|ni(^, r)| <C 1. However, inspecting on the solution (92) 
of Eq. (|87|) in absence of the Landau damping effect, 


we find that this condition may be satisfied when either 
U is small or the denominator of the expression on the 
right-hand side becomes larger. Also, for real values of 

W and U (> 0), we have TV > b' 2 when TV is indepen¬ 
dent on time r. However, when TV = TV(a, r), we must 
—-2 

have TV (a, r) > 6/ 2 Vr for Ai and A 2 to be real. Here, 

we consider the case in which TV 2 (a, r) 6/ 2 Vr. In this 
case, though some generalization may be lost, the basic 
qualitative features [21] of the solution of Eq. (109) will 
remain the same. Under this approximation the solution 


of Eq. (109) is given by a quadratic equation in TV which 

( 112 ) 


yields two solutions (for £ = ±1) 

N{a 1 r) = (3b 2 /c) 1 / 4 N 0 (X 3 /X 4 ), 
where TV = TVo at r = 0 and 

A 3 = 3 5 / 4 (6 2 /c) 1/4 TV 0 + C |"(36 2 /c) 1/2 


x (3TVo =F nbi/c) 2 ± 2\/2nbiA' clN^t/c , (113) 

A 4 = \J 362 /c ^ 6 TVo =F Znbi/cj — \f2A! oTV/Jt, (114) 

r°° r°° d dz' 

A' = P / / sechz-— (sechz') - dz. (115) 

7-OO 7- 00 oz' z-z' 
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We mention that the upper and lower signs in the ex¬ 
pressions for Ai, A 3 and A 4 are corresponding to ± in the 
Gardner solution (92). In our numerical investigation, we 


will, however, consider £ = 1 and the lower (upper) sign 
for the existence of solitons in different regions as in Fig. 
4(a) (Figs. |4(b)| and |4(c) ) to be shown shortly. Clearly, a 
decrement of the wave amplitude of the Gardner soliton 
is seen to occur. 

On the other hand, in absence of the Landau damping 
effects (i.e., when a = 0), the KdV and mKdV equations 
[Eqs. (20) and (90)] have the following traveling wave 
solutions 


= N sech 2 [(£ - Ur) /W ], 


(116) 


(z - Ur) /w] , (117) 


n = N sech 


where N = 3 U/b and N = yd>U/b' 2 are the constant 

amplitudes, W = \Jl2c/Nb = \f\cfU and W = y 1/U 

are the constant widths, and U and U are the constant 
phase speeds (normalized by c s ) of the KdV and mKdV 
solitons. 

Following the same procedure as above and assuming 
that is a small parameter and 1 ~ 6 , b 2 ~ c a 
e, analytic solutions of Eqs. (20) and (|90|) can also be 


\Z-IJ Ndrj/W +o(o), (118) 


obtained. Thus, for the KdV equation ( |20[ ) we have the 
following solution [ 2 T 

n = N(a, r) sech 2 
where 

N( a ,r) = N 0 [l + (r/r 0 )]~ 2 , (119) 

with N = Nq at r = 0 and tq being given by 


-1 a IbN 0 
T ° ~ 4 


a?-/:/; 


— 00 J —00 


sech 2 z d 
z — z' dz' 


(sech 2 /) dzdz'. 

( 120 ) 

Also, for the mKdV equation (j90l) , a substitution of b\ = 


0 in Eqs. (109)-(111) gives rise the following solution 


n = N(a, r) sech 
where 


£ -«s.r 


1-2 


N(a,r) = N 0 [1 - (r/r)] 
with N = No at r = 0 and r being given by 
aN 0 


o(a), 

( 121 ) 

( 122 ) 


r 1 = 


3\/6 V 2 


= P / [ ^^^(sech z')dzdz'. 

7_oo 7.00 z-z'dz' K ’ 

(123) 


We numerically investigate the time variations of the 
wave amplitudes of the KdV, mKdV and Gardner soli¬ 
tons given by Eqs. (118), (121) and (112) respectively. 
The results are shown in Fig. [4] for different parameter 
regimes as applicable for the existence of KdV, mKdV 
and Gardner solitons. We note that when the values of 
a and [i fall in the common regions of i?i, R 2 > 0 (for 
which 6 , bi and b 2 > 0 ), but well below the critical val¬ 
ues as in Fig. [lj the KdV equation suffices to describe 
the nonli near e volution of DA solitons. This is evident 
from Fig. |4(a)| Here, it is seen that the behaviors of the 
decay rates of the KdV and Gardner solitons are almost 
similar having their values close to each other (See the 
solid and dotted lines), however, the rate becomes higher 
(in magnitude) for the mKdV solitons. Also, the KdV 
soliton points to lower wave amplitude than the mKdV 
and Gardner solitons at large r. Thus, in the parame¬ 
ter regimes with a and /i far below the critical values, 
mKdV equation may not be relevant, however, the KdV 
and Gardner equations can give rise to approximately the 
same results for the evolution of DA solitons in plasmas. 
When values of a and /r are close to the critical values 
but in the regions i?i < 0 , i ? 2 > 0 (for which b , bi < 0 
and 6 2 > 0), the KdV solitons cease to exist (because b 
must be positive for To to be real). In this case, equations 
involving higher-order nonlinearities, namely Gardner or 
mKdV equation is necessary for the evolution of DA soli¬ 
tons. Figure |4(b)| shows that the Gardner equation is 
more appropriate than the mKdV equation as it gives 
rise to solitons with lower amplitudes than the mKdV 
solitons. On the other hand, when a and p, assume val¬ 
ues very close to the critical values but still in the com¬ 
mon regions of i?i, i? 2 > 0 (for which b , b\ sa 0 and 
b 2 > 0), Fig. 4(c) shows that the decay rates of the 


KdV and mKdV solitons get higher (in magnitude) than 
that of the Gardner soliton. In this case, KdV equation 
fails, however, Gardner equation gives better (than the 
mKdV equation) results for the properties of DA soli¬ 
tons with Landau damping. It is also found that (not 
shown in the figure) as the mass ratio m increases, the 
amplitudes N(a, r), N(a, r) and N(a, r) for KdV, mKdV 
and Gardner solitons decrease. Also, as a increases, both 
N(a,r) and N(a,r) increase, but N(a,r) decreases. Fur¬ 
thermore, an enhancement of N (a, r) and a decrement of 
N(a , t) and N(a, r) are seen to occur with an increase of 
the density ratio /i. 


VII. CONCLUSION 

We have investigated the nonlinear propagation of 
small-amplitude dust-acoustic solitary waves in a dusty 
plasma with the effects of Landau damping associated 
with the finite inertial effects of hot and cold ions. By 
considering the lower and higher order perturbations, the 
evolution of dust-acoustic solitons are shown to be gov¬ 
erned by KdV, mKdV or Gardner (KdV-mKdV) equa¬ 
tion. The existence of the critical values of <r and (de- 
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noting, respectively, the temperature and the density ra¬ 
tios of hot to cold ions) in the a — // parameter plane 
for which the KdV and mKdV solitons cease to exist 
is given. It is found that the KdV solitons with Lan¬ 
dau damping exist only for values of cr and /i well below 
their critical values. Beyond the KdV limit, i.e., when 
a and p assumes values close to their critical values, the 
Gardner equation is seen to be more appropriate than 
the rnKdV equation for the evolution of small-amplitude 
dust-acoustic solitons in bi-ion dusty plasmas. It is also 
found that regardless of the ion masses equal or not, 
there exist a common value, i.e., cr, p ss 10 for which 
both the KdV and mKdV equations fail to govern the 
dust-acoustic solitons in plasmas. The properties of the 
wave phase velocity and the Landau damping rate are 
also studied with different values of cr, /r and the mass 
ratio m of cold to hot ions. It is shown that the damping 
rate is quite different for dusty plasmas with equal and 
unequal masses of positive ions. Distinctive features of 


the KdV, mKdV and Gardner solitons with a small (but 
finite) effect of the Landau damping are compared for ap¬ 
propriate parameter values of cr, /i and m. It is seen that 
beyond the KdV limit, Gardner equation always points 
to smaller wave amplitudes with the effect of Landau 
damping. The results should be useful for understanding 
the collisionless wave damping of dust-acoustic solitons 
in laboratory as well as space plasmas such as those in 
planetary rings (e.g., the F-ring of Saturn) where nega¬ 
tively charged dusts and two species of positive ions are 
the major constituents. 
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